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I. INTRODUCTION 



1»1 Background 

Although extensive ocean ambient noise studies were conducted 
dviring World War II with continuing effort since that date, few of 
either theoretical or experimental treatments have dealt with the 
directional propei*ties of such noise. There are many possible sour- 
ces of ambient noise in the ocean any of which could produce non- 
uniform dietributions. R. A. Prosch^^^ in examining likely contrib- 
utors to the total acoustic energy content of the oceans has con- 
sidered, for example, biological soxirces, seismie noise, shipping, 
lightning and sea surface waves. His conclusion that sea surface 
waves are a major contributor is in agreement with World War II 
data. For example, the "Knudsen"^^^ curves, a compilation of onmi- 
directional measurements relating ambient noise spectrum levels to 
frequency as a function of sea state, exhibit a direct sea state/ 
ambient noise interdependence and by their very nature permit the 
inference of noise anisotropy in the ocean. 

Theoretical treatments considering the ocean surface as a 
possible noise source have been published by Eckart^^^, Roberson^^\ 
Urick^^^ and Miller^^^. Eekart has presented the physical picture of 
the ocean surface resolved into sinusoidal travelling waves of the 
particular phase velocity which would be associated with free flex- 
ural waves of an elastic plate excited by random local shocks of 
short duration. The above model would result in the generation of 
sound only under the unlikely condition that the surface phase 
velocity exceeded the velocity of the bounded medium, Roberson 
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without attention to a phyolcal model of noise production considered 
both an infinite volume distribution of simple sources and the com- 
bined effect of a volume distribution of simple sources bounded by a 
surface distribution of dipole sources. Assuming an attenuation law 
proportional to frequency squared^ he showed a -6 decibel/octave de- 
pendence of noise upon frequency for the unbounded case and a non- 
uniform dependence upon frequency for the mixed ease depending upon 
the ratio of the contributions from the two sources and the fre- 
quency. In this latter case, 'no combination of surface and volume 
sources produced a constant slope In agreement with the "Knudsen" 
curves except over a limited frequency range. 

Urick specifically considered the subject of directional 
distribution of ambient noise and concluded that significant aniso- 
tropy might exist depending upon the unlcnown manner in which each 
element of surface radiates. 

Miller considered the directional distribution resulting 
from the combination of straight path propagation from surface rad- 
iators close to the receiver and convergent zone propagation from 
distant surface radiators. The acoustic radiation from the surface 
varied as cosine d , where ^ is the angle between the direction of 
radiation and a normal to the radiating surface. No bottom reflec- 
tion or volume sources were assvaned. He concluded, tentatively, 
that noise in the above case is essentially isotropic over a hemis- 
phere with a slight Increase within + 9 degrees of the horizontal 
and that the area immediately above and close to the receiver 
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contributes about the sane total noise as the contribution from the 
first convergence zone. 

If the ocean surface lea major contributor to ambient noise, 
it has been reasoned that ambient level should decrease with depth. 
The change in ambient noise levels with depth^^) has been investi- 
gated experimentally. No such change has been detected to depths of 
500 feet in tests at 25 kilocycles per second, Urick^^^ has shovm, 
however, that the change to bo expected between 5 ^ind 5OO feet depth, 
even at a higher frequency of’ 50 kilocycles per second where attenu- 
ation is greater, is only 5 decibels theoretically and thus a smaller 
change would have been obscured by measurement fluctuations at the 
lower frequencies, Teer^®) reported the results of teats at 19 
kilocycles per second vfith a small line array oriented horizontally 
and vertically in deep water. The results, which seemed to indicate 
that ambient noise is isotropic, are not considered conclusive be- 
cause of limitations In equipment and the environmental conditions 
experienced during the tests. 

In summary, no experimental data exist to describe the di- 
rectional distribution of ocean ambient noise. Theoretical treat- 
ments have indicated either laotropie or anisotropic distributions 
may exist depending upon initial assumptions as to the nature of the 
surface radiator and the uncertain effect of ocean sound channels. 

1 ,2 Research Objectives 

It is the purpose of this study to investigate anisotropy 
of deep ocean ambient noise in the frequency range of 750 "to 1500 
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cycles per second. Further specific objectires will bet 

a. To isolate rertical and azimuthal dependencies, 

b. To establish the identity of the sources contributing 
to non-uniform distributions. 

e« To infer the radiation patterns of such sources. 
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II, AMBIENT NOISE EQUIPMENT 
2,1 Dlaeusfllon 

Thera are two major reaeone which hare prevented more wide- 
spread investigation of the directional properties of ambient noise. 
The first concerns the limitation of a single directional hydrophone 
which, to monitor several directions in space suceessively, reqvdres 
the assixmptlon of a stationary proeess. The above difficulty may be 
overcome by the employment of more than one directional hydrophone 
or by forming several directional beams from an array of fixed non- 
directional hydrophone elements. Either solution generates a second 
problem, namely the difficulty of ensuring unifora gain in the sev- 
eral beams. Both limitations have been surmounted by a system con- 
ceived by Anderson^^^ and Implemented by the staff of the Marine 
Physical Laboratory of the Seripps Institution of Oceanography, 
University of California, La Jolla, A 52-element array has been 
constructed in the shape of a great stellated ioosohedron as shown 
in Fig, 1, A summary of the array design parameters is contained in 
paragraph 2,2 following. The 52 hydrophone outputs are combined to 
fora 52 beams with axes coinciding with lines passing through the 
array center and the Individual hydrophones. All 52 hydrophones are 
used in forming each beam. 

The problem of uniform sensitivity in each beam is minimized 
by employing clipped processing. At the output of each hydrophone 
channel prior to beam formation, all amplitude information is dis- 
carded by infinite clipping and phase information only retained. 

The clipped "on-off" form of the channel outputs permits digital 



5 




FIGURE I. THE AMBIENT NOISE ARRAY 
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shift register techniques to be employed in the forming of the 52 
beams with an attendant saying in space and equipment complexity 
oyer conyentional analog beam forming techniques. The system is 
thus independent of indiyidual channel gain differences from each 
hydrophone thro\igh the beam forming process. The problem of ampli- 
tude differences appears only in compensating for detector gain dif- 
ferences in the process of rectifying each beam output. The clipping 
and beam forming processes utilized in this system hays been re- 
ported by Anderson^ and Rudniek^^^^. The restoration of the ana- 
log form of a signal from clipped hydrophone oxitputs is considered 
briefly in paragraph 2,5, 

2.2 Array Configuration Deyelopment 

In arriylng at the present array geometry of a great stel- 
lated icosohedron, Anderson^^^ considered three factorsi 

a. The number of elements. 

b. Element configuration. 

0. Element spacing. 

Since presumably small deviations from an isotropic field were to be 
measured, a symmetrical distribution was suggested permitting nearly 
identical directly ity patterns of the several beams to be formed by 
electrical steering. A distribution of elements in depth rather 
than a spherical shell would permit greater interelement spacings 
for a given array diameter and thus increased array efficiency 
through minimizing Interelement correlations. The number of elements 
was limited by considerations of electronic complexity. Thirty-two 
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appeared a practical nmber and (a) and (b) above were satisfied by 
the great stellated icosohedron. Optimun design of intoroleaent 
spacing was detomined by calculations of array efficiency defined 
as the ratio of output to input signal to noise ratios. Signal in 
this case vjus the deviation from isotropy which it was desired to 
measure in a given direction. An interelement spacing of 0.6 wave 
lengths at the 1 kilocycle per second nominal center frequency vreis 
chosen. 



2.3 General System Description 

An overall system block diagram is shoim in Fig. 2. For 
simplicity only 1 of 52 identical receiving channels is represented 
completely. Reading from the left, an individual hydrophone output 
is first aEq)llfied and then employed to frequency modulate a tele- 
metering transmitter. The 52 transmitter carriers are distributed 
from 10 through 16.2 megacycles per second in 200 kHocyoles per 
second increments. An additional 1 megacycle per second amplitude 
modulated tranoaitter relays hydrophone array azimuthal orientation 
in accordance with the directional indication of a magnetic compass. 
The 52 hydrophones, preamplifiers, telemetering transmitters and 
compass system comprise the undenrator portion of the system. 

The transmitter outputs are coupled to the ihlpbome re- 
ceiving portion by coaxial cable. A I50O foot length was available. 
The 5/8 inch diameter, doubly annored cable seirves also as the 
supporting member of the 700 pound receiving array. Cable attenua- 
tions of 1.5 and 1.8 decibels per 100 feet at 10 and I6 megacycles 
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FIGURE 2. AMBIENT NOISE EQUIPMENT BLOCK DIAGRAM 
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per second respectively necessitated the insertion of a broadband 10 
to 1^,2 cierjacycle per second preeanplifler prior to deaodulation of 
the signal by the J2 ^latching frequency modulated receivers. At this 
stage It is possible to vlev; the amplitude and spacing of the J2 in- 
dividual carriers as irell as the frequency deviation imposed by each 
signal by means of a panoramic receiver and cathode ray tube indica- 
tor. 

Following demodulation in the FM receivers, the audio signal 
is further amplified and then infinitely clipped. All amplitude in- 
formation ie discarded and only polarity information retained. The 
clipped ttgnal then enters a digital shift register ten bit delay 
line which provides suitable signal delay Increments such that sub- 
sequent ecabination of the signals in a mixing matrix yields J2 
beams symmetrically spaced over 4 7tS2» The J2 beam outputs are then 
individually detected and Integrated by a simple low pass RC filter 
of one second time constant. Following detection and integration, 
the p2 beam outputs are sequentially sampled at a 50 cycle per second 
rate by a mercury jet scanning switch for presentation on a cathode 
ray oocillosoope. Sixteen millimeter scope photographs preserve the 
data for analysis. 

All receiving equipment is housed in a portable Army Signal 
Corps communications hut. The entire 7|000 pound installation may 
be moved quickly by crane from Its laboratory position to the fantall 
of PAOLIKA T, a Scripps Institution of Oceanography research vessel. 
The electronic hut, array, associated winches and boom end gasoline 
powered 15 KW generator satisfying the I 15 volt, 50 ampere demand 
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are Bhown in Figs. 5 

2,4 Hydrophone and Transmitter Assembly 

The hydrophones are encapsulated barlun titanate cylinders 
with a flat response averaging 90 decibels below 1 volt//^ bar 
throu^out the range of Interest, 750-1500 cycles per second. The 
cylindrical sensitive element is 1-1/2 inches high, 1-?/l6 inches in 
diameter and has a wall thickness of l/8 inch. Individxial hydro- 
phone directivity is discussed in paragraph 5»5* Although the par- 
ticular unite employed were constructed at the Marine Physical Lab- 
oratory, an identical unit, the REMAOO Model R-15O is ccomeroially 
available. They are rated at a working pressure in excess of 16OO 
pounds per square inch. The hydrophone fits into the end of a 
16-1/2 inch long, 5/52 inch thick stainless steel pipe which forms 
the watertight housing for preamplifier and transmitter. All h3rdro- 
phone units are series connected, 50 VDC reqvdred by each unit or 
1600 VDC total. A relay feature is included to remove a particular 
hydrophone unit should a vacuum tube filament open. Were this fea- 
ture not included the entire array would become disabled by the 
failure of a single tube. On the other hand, a short circuit any- 
where in the series string disables all units beyond the point of 
shorting, A partially disassembled unit is shown in Pig. 5* 

2 <5 Clipped Processing 

Because digital multibeam steering (DIMUS) is a relatively 
now technique, the following simplified illustration of the restora- 
tion of the analog form of a signal by the addition of several 
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FIGURE 3. AMBIENT NOISE EQUIPMENT 
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FIGURE 4 AMBIENT NOISE ARRAY WINCH 





FIGURE 5. HYDROPHONE / TRANSMITTER ASSEMBLY 
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eamples of the signal which have been clipped is included. Referring 
to Fig, 6a, the time varying analog audio output of an individual 
hydrophone channel is shovm. Some band limiting haa occurred in the 
hydrophone preamplifier bandpass filters. Further high-pass filter- 
ing by the simple network of Pig, 6b removes any remaining fluctua- 
tion below the frequency of Interest resulting in the waveform of 
Fig, 6c, The waveform of 6c is transformed into one of two states 
as shown in Fig, 6d by the clipper amplifiers. Figure 7a repeats 
Fig, 6d, In addition Figs, 7b, ->c, d and e represent the outputs of 
four other channels. The example has been so chosen as to illustrate 
the array gain of a five hydrophone system in detecting a directional 
signal mixed with omnidirectional noise. An input signal -to-nolse 
ratio at the individual hydrophone of xmity has been assumed. The 
direct addition of the five output voltages produces a randomly 
fluctuating voltage as shown in Fig, 7f« The absence of a prcmounced 
peak is to be expected since the noise at each hydrophone due to the 
isotropic background is uncorrelated while correlation between volt- 
ages generated by the directional noise will be visible only when 
the proper electrical delays corresponding to the geometrical delays 
between hydrophones are introduced by the shift registers. The 
voltages of Pigs, 7a through c are sampled at the shift register in- 
puts at a 5 kilocycle per second rate. This timing satisfies the 
Nyquist criterion of a minimum sampling rate equal to twice the 
highest frequency of the spectrum of interest. Figures 7g through 
7k repeat Figs, 7a through 7e with the proper delay introduced by 
the shift registers to correct for hydrophone geometry. The 
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FIGURE 7 THE ADDITION OF CLIPPED SIGNALS 
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addition of Figs. 7g through 7k shown in Fig. 71 again produces a 
randomly fluctuating voltage with the exception that now a peak ap- 
pears due to the directional noise. The reproduction of the original 
signal improves with the number of channels employed. In aural test 
with the 52 element system it is not possible to distingitiBh between 
an undipped audio signal input as monitored at the output of a sin- 
gle hydrophone and the reconstituted input as heard at the output of 
a beam formed from the clipped signals of the 52 elements. 

2.6 Array Beam Patterns 

The sunmed voltage output of all hydrophones may be written 1 

%ut^Z e [cos(^- + + 

/-/ ^ L 

where u •= angular frequency, radians per second 
“ hydrophone coordinates 
Og, m electrically steered direction 

“ angular deviations from steered direction 
Due to hydrophone symmetry, a particular beam will exhibit conical 
symmetry about the beam axis. The ioosohedron has two axes of sym- 
metry and, accordingly, two different beam patterns. In a spherical 
coordinate system with the array center as origin, every element can 
be pictured as positioned on one of two concentric shells. There are 
20 elements located on the outer shell and 12 elements on the inner 
shell. The two axes of symmetry are formed by lines through the 
array center and elements on the outer and inner shells. Since each 
array element has a diametrically located twin, the expression for 







. 




19 



power output with = 0 and = 0 becaaea; 

* ^2 r D r 2 

G {0q) = J cos cos^z^. |cos - cos( ) j )| 5* (^y ) dou 

where io a v;eighting function, the relative frequency recponae 

of the system and a>^, a>^ are the system bandpass limits. Computed 
patterns normalized to unit on-axis response are shovm in Fig. 8. 

The continuous curve represents the solution to the preceding equa- 
tion for an outer beam (a beam whose axis passes through an outer 
element in the array). ^ 

The use of the preceding equation as written tacitly assumes 

that the exact time delay specified by the expression COS<^ 

c 

is arailable. In praoties tha exact delays ere not available owing 
to the llnitations of the digital shift register which at its tapped 
outputs provides delays equal to integer multiples of the clock gen- 
erator period. The dashed curve in Fig* 8 is the outer beam re- 
sponse computed on the basis of the aetml shift register delays. 

The dotted curve ie the corresponding response for an inner beam. 

The ecmblned effect of beam differences due to element configura- 
tions and steering delays results in an inner to outer on-axis beam 
response ratio of 0.8^^, a factor which must be considered in the 
data reduction process. 

2,7 Reduction of External Noise Interference 

The preliminaiy phases of at sea measuremente confirmed the 
existance of two eoturces of interfering noise which had been antici- 
pated, They were; 
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a. Ship nachinery noise tramsmltted through the hull struc- 
ture into the ocean. 

b. Self-noise generated by relative ooveraent of the array 
and ocean water. 

Another troublesome source, external shipping interference, could 
only be controlled by judicious choice of operating area free from 
shipping lanes or by testing during night tine hours when shipping 
traffic was at a mininun. A fourth source generated by the rolling 
motion of the testing ship is discussed in paragraph 

Own ship machinery noise was eliminated by two means. First, 
all normal ship machinery was secured during the testing period. 
Secondly, interference from the special 15 KW gasoline powered gen- 
erator serving as the power supply for the ambient noise equipment 
was reduced to a negligible value by isolating the generator from 
the ship deck by autcanobile tire casings. The effectiveness of the 
isolation method was demonstrated by comparative evaluation with a 
special 5 KW sound isolated unit developed at the U, S» Navy Elec- 
tronics Laboratory for similar noise teste. This latter source was 
completely suspended from each of its four comers in a pipe support 
frame by airplane shock cord. Although its electrical capacity was 
not adequate for full time employment, alternate operation with 
first the limited capacity 5 KW unit and then the larger I5 KW 
source revealed no differences in ambient noise levels as measured 
at the array suspended in the ocean. 

Elimination of the second interfering source, self-noise 
generated by relative array/water motion, was accomplished by two 



nepjis: (1) the suapenolon of the array fron a eyetem of buoys and 

floats In order to isolate the array from ocean surface and ship mo- 
tion and (2) the employment of a parachirte sea anchor to eliminate 
ship drift caused by s^urface winds and attendant array/water rela- 
tive motion. The two techniques are shown achematically in Fig. 9« 
The suspension was so designed that the spar buoy would float in the 
position shown independent of array depth. The floats compensated 
for cable weight as the array was positioned at 100, ^00 or 1000 foct 
depths. The spar buoy acts effectively as a spring with constant, 
k , equal to the buoyancy force per length submerged. For the A inch 
diameter buoy in question A" = 5 lbs/ sec^ which for an array dis- 
placement of 5^0 pounds results in a natural period of oscillation of 

T - » 60 seconds 

V k 

a value 10 times greater than the longest period of the ocean waves 
encoiaitered during the experiment. Particular attention was devoted 
to buoy swivels and connecting links in order to prevent noise gen- 
eration by any metal-to-metal contact. 

The sea anchor employed was a siurplua government parachute 
fastened in the manner shown by manila hawser and equipped with a 
manila tripping line for collapsing purposes. Even in very light 
wind of 7 knots it was effective In maintaining ship's head into the 
seas, reducing ship roll end completely eliminating drift. Under 
certain not Infrequent conditions stratified ocean currents were ex- 
perienced where surface currents and wind were in a direction differ- 
ent from the ocean current at the array depth. The net effect was the 
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SEA TEST CONFIGURATION 



towing of the array by the test ship and parachute and invariably 
resulted In test cancellation. 



III. SYSTEM MEASUREMENT ACCURACY 



5*1 Discussion 

The ideal manner in vhieh to specify system, measurement ac- 
curacy vjould be in terras of beam response to a calibrated source cap- 
able of producing an output spectrum covering the ambient noise 
equipment bandwidth. Up to the present the above procedure has not 
been possible for a variety of reasons. The chief limitation arises 
from the low frequency of operation of the system. Neither the nec- 
essary low frequency broadband source nor a qidet, deep ^mter area 
with adequate facilities for maintaining array/source geometry has 
been available for test. Although it has not been possible to cali- 
brate the system as a single unit by one series of measurements, 
calibration of all component parts of the system, individually, has 
been possible and is considered an adequate substitution for the 
single measurement procedure. In paragraphs 5*2 throiigh J.7 follow- 
ing, the possible error contributions of all component parts of the 
system are examined. Detector calibration has been established as 
the major source of error In the system and will accordingly be oon- 
sidered first. It is the objective of the remaining paragraphs to 
show that all contributions from other components of the system are 
negligible in contrast to detector calibration error. 

3.2 Detector Calibration 

The use of clipped processing considerably ccmplicated de- 
tector calibration procedures. Since the output amplitude of a 
given beam is a fxmction of the relative phase correlation between 
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Individual hydrophones and Is Independent of individual hydrophone 
output amplitude, the input/output charaeteristio of a given beam 
foming/detector combination could not bo established by Injecting a 
signal of controlled amplitude Into each clipper amplifier and mea- 
suring detector output voltage. 

The problem Is fuller complicated by the nature of the shift 
register matrix output. A representative waveform in shown in Fig. 
10 a lllustiratlng the summed voltage oirtput from 52 shift registers 
with random noise input. The desired information is contained in 
the fluotvtating AC. The mdesired 5 kilocycle per second clock 
pulses, arising from the particular form of shift register circuitry 
employed, eliminated the possibility of mploying full wave rectifi- 
cation and necessitated use of less efficient half-wave rectifiers. 

A quantitative analysis of the system output Including the effects 
of clipping is contained in Appendix I. 

Although 52 random noise generators whoso relative conrela- 
tions ootdd be varied at will would permit detector calibration, no 
circuitry to fulfill this requirement was obvious. Instead a per- 
fectly correlated AO signal was employed. A 2.5 kilocyele per secoid 
tost voltage of two phases, 180 degrees apart, was derived by fre- 
quency divider from and synchronised with the 5 kilooycle per second 
clock generator. Application of one or the other of the two phases 
to a partlcxilar clipper amplifier input would assure that the corre- 
sponding shift register outputs were identical, either off or on. 

By adjusting the number of clipper amplifier Inputs of each phase, 
the matrix output for any beam could be varied from 0 to 4 16 units. 
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OSCILLOSCOPE TEST PATTERN 
10(c) 



FIGURE 10 DETECTOR OUTPUT PATTERN 



Detector bias was adjusted so that vrith 17 Inputs of one phase and 
15 of the opposite phase, the detector output amplitude was 1/2 the 
difference in clock pulse amplitudes as viewed by oscilloscope and as 
shown in Fig* 10b» Transient spikes as indicated in the figxire and 
detector inefficiency prevented the achieving of exactly zero output 
voltage for the condition of inputs divided equally as to phase* 

Time variations in individual detector characteristics were 
compensated for by trimming each detector bias to permit alignment 
with a predetermined test pattern as shown in Pig. 10c, This was 
created by steering the test si^al to produce maximum response c»i 
beams 1 and J2, 

By progressively changing the inputs from half of each phase 
(minimum deteotor output) to all of one phase (maxiram detector out- 
put) complete detector characteristics could be traced, Meuiy trac- 
ings for each detector were required to average out the slight dif- 
ferences in individual register outputs emphasized by the artificial 
nature of the single frequency input but averaged automatically In 
practice by the random nature of the noise. The applicable porti<m 
of a typical oiurve is shown in Fig. 11. After one set of character- 
istios had been determined as described, a duplicate sot was pro- 
duced, In no case did the second set differ by more than 1 millivolt 
from the first set. This 1 millivolt error is considered the maxi- 
mum system error and amounts to an error in the measurement of the 
output of any beam of 0,2 decibels. 
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FIGURE II TYPICAL DETECTOR CHARACTERISTICS 



50 



5» 5 The Array 

Individual array olenents compricing the BaTiO^ receiving 
unit and the FM transmitter canister housing were calibrated by the 
U. S. Navy Electronics Laboratory calibration facility at Sweetwater 
Lake, San Diego, California, Measured horizontal and vertical re- 
ceiving patterns were completely circular. Since the v;avelength at 
the nominal array center frequency of 1 kilocycle per second is 5 
feet, the patterns are understandably uniform indicating that the 
individual hydrophones are equivalent to a point hydrophone. The 
cylindrical shape of the hydrophone and canister does not alter the 
isotropic response. 

Although the response of the assembled array of individual 
elements could not be measured, acoustic theory and V/orld \’ar II 
measurements^^ of sonar deme structures demonstrated that support 
members of cross-section equal to or less than I/15 of a wavelength 
have no effect upon the beam patterns of arrays transmitting through 
such structures. Since the ambient noise array support members are 
1 inch steel angle Iron approximately I/8 inch thick, the above de- 
sign criterion is more than adequately satisfied. It is concluded, 
therefore, that the array may be considered fojrmed of exactly 
identical annidircctional receiving olamonts unaffected acoustically 
by any supporting structure. Accordingly, it is considered that no 
errors are Introduced by the array, 

5*4 The RF System 

The RF system was examined as to circuit noise, dynamic 
range and the introduction of undesired phase shifts. 
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With an appropriate capacitor aubgtituted for the hydrophone 
cryatal, electronic circuit noise levels averaged 20 decibel e below 
sea state zero noise levels. 

The question of dynamic range primarily concerned channel 
cross-modulation. Modulation level* wore adjusted *ueh that maximum 
noiee levels as measured at eea produced no more than a + 50 Icllo- 
cyole per second deviation of the appropriate channel carrier. This 
adjustment permitted a 100 kilocycle per second eeparatlon between 
channels and ellodnated all possibility of cross-talk between chan- 
nels. 

Signal phase shifts unintentionally Introduced by receiving 
channel components are equivalent to electrical time delays. Only 
differences between channel phaee shifts are significant since equal 
phase shifts in each channel are equivalent to the application of a 
fixed time delay to all channels which would not affect beam output. 
Channel phase shift differences will not seriously affect beam out- 
put until they approximately equal the sample period of 0.2 milli- 
seconds or 72 degrees. A large relative phase shift would cause the 
signal to appear at the incorrect shift register tap and would dis- 
tort the beam ppttem. To measure the actual shift introduced by 
the system from hydrophone to shift register Input, conventional 
Llssajous figure osoillosoope techniques were employed. An audio 
oscillator output was coupled by suitable attenuators to the osolllo- 
scope horizontal deflection plates. A parallel oscillator output 
was capacltatively coupled to the hydrophones in air by a special 
cup designed to fit concentrically over the array hydrophone 
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elements. The clipper amplifier output of the particular hydrophone 
channel undergoing measurement was fed to the vertical deflection 
plates of the oscilloocope. The oscillator frequency was adjusted 
for a straight line or zero phase shift. For the channels the 
frequency of zero phase shift ;7as fomd to vary between 4_ 50 cycles 
per second centered at 1,01 kilocycles per second. A second measure- 
ment of the phase shift characteristic of an individual channel 
determined by observing the frequency corresponding to the easily 
read Liosajous points of 180, 90 and zero degrees permitted conver- 
sion of the initial zero phase shift frequency measurement to phase 
shift at zero frequency. Front this latter auxiliary measurement it 
was established that 0.55 degrees shift occurred per cycle or that 
the greatest measured shift was + 10 degrees from 1.01 kilocycles 
per second, a value well within the sampling period. 

It is concluded, therefore, that no significant errors tirere 
introduced by the RF system, 

5.5 Beam Forming System 

Errors in beam output introduced by the beam forming system 
may occur due to four possible causesJ 

a. Inadequate clipper amplifier dynamic range. 

b. Nonimiformity of clipper amplifier input band-limiting 

filters. 

c. Amplitude variation of shift register outputs, 

d. Nonuntfornlty of matrix summing resistors. 
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Clipper amplifier dynamic range, defined as that range where 
clipped output is independent of input amplitude, was measured to be 
70 decibels, a value well in excess of the range of ambient noise 
levels in the ocean* No error is attributed to this source. Ander- 
son(^5) has shown that errors in output of Individual elements com- 
prising an array degrade the array signal -to-noise output ratio in 
accordance with the relationJ 

S/N = /7 { I - ^7-2) 

where n is the number of array elements and cr is the standard devb- 
tlon of the individual element errors. The standard deviation for 
items (b), (c) and (d) above when measured Indlvldiially and then 
combined was less than 10 per cent. Substitution of this value into 
the above expression yields a 0,05 decibels loss, a negligible con- 
tribution to overall system error from beam forming ccxnponents. 

3*6 Film Recording 

As noted in paragraph 2.5* Individual beam outputs were 
photographically recorded continuously throughout each run by l6mm 
camera. Frame spacing was geared down to 5»5 frames per second. The 
developed film was viewed for reading on a Craig l6mm Projecto Edi- 
tor. The displayed frame size was as shown in Fig, 12 where a ^ 
millivolt calibration pedestal, the base line and 52 beam outputs 
are Indicated, An accurate, fine line, overlay grid of h millivolts 
per division was drawn to large scale and then photographically re- 
duced to the size of Fig, 12. The reading procedure first Involved 
adjusting the projected frame size such that the calibration pulse 
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FIGURE 12. FILM RECORD AND OVERLAY MEASURING GRID 
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height corresponded exactly to AO nillivolta on the viewing grid. 

The value of each beam output could then be read quite ainply to the 
closest millivolt. Assuming that individual component errors add 
randomly, the error introduced in the film reading process is con- 
sidered small in comparison to the detector calibration error of 
paragraph 5*2. 

3«7 System Linearity 

A consequence of digital array phasing by the DIMUS technique 
and symmetry in the array beams, is that the sum of the output powers 
of the 52 array beams is a constant. Thus, an Inerdaee in level in a 
given beam will result in a corresponding decrease in level in the 
remaining beam outputs, A very strong signal, or its equivalent, 
good correlation in a given direction, can result in beam capture 
with reshltant suppression of the remaining beam outputs and attend- 
ant loss of information as to the relative field intensities viewed 
by several beams. It is necessary, therefore, to determine that 
maximum output level on a given beam which does not introduce dis- 
tortion in the outputs of the remaining beams. Toward that end it 
is assmed that a signal-toWnoise ratio at the input of an individual 
hydrophone of 1 to 10 represents a sufficiently small signal such as 
to preclude the po3sibllit3'- of signal capture. The output power gen- 
erated by such a signal-to-noise ratio has been derived in Appendix 
II and the results compared to maximum outputs observed during mea- 
surements at sea. Since no measured beam outputs exceeded those 
computed, it is concluded that the system behaved in a linear manner 



throughout the experiment 
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IV. TEST PROCEDURE 
4,1 Narrative 

The on-station test procedures were straightforward barring 
the occasional equipment component failure. Two hours prior to 
ccmmencement of test the equipment was energized in order to permit 
a stable temperature to be reached. The detector circuits, for ex- 
ample, proved to be extranely sensitive to temperature changes 
greater than 1 degree and required thermocouple control of the two 

electronic equipment hut cooling fans. Upon reaching the test area 

- 

the parachxite sea anchor was rigged and the ship permitted to reach 
an equilibrium condition headed into the seas. The array was then 
swung clear of the ship and lowered to within approximately 80 feet 
of the desired operating depth. The spar buoy and float system was 
then damped to the array cable and the system lowered until the 
weight of the array was supported by the floats. The array was po- 
sitioned generally between 100 and 150 feet from the ship. 

Radar bearings and ranges to any oontacts, anemometer read- 
ings and fathometer soundings were recorded prior to turning off all 
ship power with the exception of the ambient noise equipment genera- 
tor, The detectors were balanced to the calibration signal as des- 
cribed in paragraph 5*2 and an approximate 4 minute film recording 
made of backgroxmd noise. Array ccmipass heading, ship's head, swell 
direction, bearing to the array, oscilloscope and calibration pedes- 
tal and base line settings were recorded. 

Normally 5 recordings over an hoiu* interval were made at 
each depth. At the conclusion of a given sequence of recordings at 
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one depth, ship's power equipment was energized, the array partially 
raised to permit detaching the float system and then lowered to the 
new depth. The procedure as described above was then repeated. Ap- 
proximately 6 hours were reqtdred to complete measurements at one 
station of 5 array depths. 

4.2 Test Environment 

Tests at sea were planned with the objective of maintaining 
a constant environment. Since ocean surface effects as noted In 
Section I were a likely source ot non-uniform ambient noise distri- 
butions, an effort was made to control other possible variables such 
as the effects of land masses, surf, shallow water and varying types 
of ocean bottom. Although It would have been desirable with the 
above objective In view to c(»iduet all tests In the truly deep ocean 
beyond the continental shelf, only one such operation was undertaken 
owing to the limitations In ship speed which dictated a 24 hour one- 
way travel time to such areas. Only one such extended cruise was 
undertaken and accounted for all sea state three, 2000 fathom ocean 
depth data to be reported. The sea state one plus data were obtain- 
ed In 1000 fa them water In an area which shoaled to 500 fathoma In 
all directions within approximately 10 miles. The sea state one- 
half data were obtained In a 1000 fathom area similar to that 
studied vnder sea state one plus. 

Sea state classifications are at best an approximation even 
when derived from a stable platform. Since the difficulty of main- 
taining an objective viewpoint from a rolling and pitching 75 foot 
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tort platfom was appreciated, the initially recorded sea state 
values of two and fotir were reduced to one plus and three respec- 
tively following examination of motion pictures taken at the time of 
recording. All observations wore conducted in moderate swell con- 
ditions, The environmental conditions are siaamariaed as follows: 



Sea State 


Ocean Depth Ajrray Depth Type Bottom 


5 


2000 fathoms 150, Gray mud 

1000 feet 


1 + 


1000 fathoms *- I50 feet Gray mud 


1/2 


1000 fathoms 1000 feet Gray mud 



7. MEASURED DATA ANALYSIS AND RESULTS 



Noise Sovurce Origin 

Nan-iBotropic noise distributions, functions of both azimuth 
and elevation angle, ware measiu*ed during the tests described in 
paragraph 4,2. Observations of the character of the noise fluctua- 
tions in the upper beams of the array, namely the rertieal upward 
beam directed normal to the ocean surface and the next family of 5 
beams at degrees from the vertical, indicate that one of the sev- 
eral noise generating mechanisms is the white cape of breaking 
waves. Near the ocean surface individual beam outputs vary consid- 
erably from one another in a random manner. Momentary increases in 
level in a given beam were probably duo to the arrival of energy 
from an individual or small group of white caps or wavelets. As 
array depth increased, boam-to-beam differences decreased. The 
logical explanation lies in the Increased area of ocean surface sub- 
tended by the beams and the likllhood that the beams of a particular 
elevation angle view more nearly the same distribution of white caps. 
The above ooneluslon was substantiated during tests designed to re- 
duce external noise interference described in paragraph 2,6, 

Despite efforts to sound isolate the 15 KW power generating 
equipment during early phases of the experiment. Interference from 
the test ship invariably was observed in that beam most nearly di- 
rected at the ship. When installation of special sound isolated 
generating equipment failed to reduce the interference, further ex- 
amination revealed that the test ship was vmwlttlngly generating the 
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oquivalont of white caps. The heavy handling boom, rigged to port 
for all tests oreated a pronounced list. The list, coupled vriLth a 1 
to 2 foot freeboard aaidships, resulted in the deck edge dipping be- 
low the waves as the ship rolled even in the lowest sea states. A 
considerable voluae of water entered through deck edge scuppers. On 
the return roll this water would pour from the scuppers on to the 
ocean surface like water from a pitcher. The process repeated 
with each roll. The intei’ference was substantially reduced and the 
source identified thereby by streamlining the entire portside, 

■a 

rounding over projections and covering scuppers with a 40 foot long, 
4 foot wide strip of canvas. 

White caps are probably the principal contributor to back- 
grovmd noise. However, the ocean surface by some vmexplained mech- 
anism still contributes to the noise field at sea states above sero 
even when white caps are not present. This fact is eubstantlated by 
the measxjired distributions reported in paragraph 6.5. 

3,2 Data Analysis 

The film data wae first reduced to tabvilar form as described 
in paragraph 5,6, Three equally spaced 4 minute samples recorded 
within a 1 hour Interval and with sin observation every 14 seconds 
within each sample resulted. The 5 ssunples for each array position 
were examined for fluctuation* within each sample and between sam- 
ples. Fluctuations within samples varied from lees than + 0.4 deci- 
bels for beams directed in the lower hemisphere to + 1,0 decibels 
for beams directed upward. The standard deviation of sample 



averages was smaller generally than the standard deviation of obser- 
vations within a sample by the factor y'^ where n is the 

ntanber of observations ocmiprlsing a sample. Accordingly three sam- 
ples at each depth were averaged and treated as one, 

^.3 Azlmtrthal Dependence 

The data were first analyzed for azimuthal dependence by ex- 
amining the outputs of families of beams of a particular elevation 
angle. In general, some azimuthal dependence was apparent princi- 
pally in the outputs of those array beams directed in the upper 
hemisphere. The distributions normally were elliptical in shape, 

A representative azimuthal distribution is ahm-m in Fig. 15* It is 
considered significant that the maximum beam output for a given ele- 
vation angle normally aligned Itself with the wave fronts of the 
swells while masrtmm outputs usually were perpendicular to the 
swells. This fact suggests a non-uniform radiation pattern of the 
white cap or wavelet source; transmission in directions normal to 
the ^vave front is impeded by the troughs, with lateral transmission 
along the ridges being less affected, 

5*4 Vertical Dependence 

The data were examined for non-uniformity in the vertical 
plane by averaging the outputs of all booms of a given elevation 
angle, thus eliminating azimuthal dependence. The results are shown 
in Fig, 14 which presents the measiured field as a fiaiction of eleva- 
tion angle for the conditions listed in paragraph 4,2, Sea state 
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effects are apparent v.dth maximum to mininim: differences for a given 
condition varying from less than 1 decibel for sea state one-half to 
6 decibels for the highest sea state. Less significance can be at- 
tributed to differences in anisotropy magnitude due the variable of 
array depth because of mcertalnty as to the constancy of sea state 
during the recording period. Approximately 6 hours were required to 
record data at the three depths during an early morning period of 
variable vfinds. It is not unlikely that the recorded sea state of 
three could in fact have fluctiiated between sea state two plus and 
three plus. 



n, TRIE FIELD RESTORATI<^ 



6,^ True Flold Indeterminacy 

In the preceding paragraphs the analysis has been confined 
to the measured fields. It is necessary to distinguish between the 
measured distributions and the true noise distributions. The output 
from a beam whose axis is oriented in a fixed direction is not de~ 
termlned solely by the true field intensity in the oriented direc- 
tion. The beam output is proportional to the weighted mean of the 
intensities in all dlreetions. The weighting factor is the beam 
pattern. Mathematically expressed, the output of a beam becomess 



where / ■ beam output or measured field. 

G = the beam response, a function of the angular difference 
between axis orientation, dg , ■, and a direction of 

interest ^ 

X » the true field intensity where, for simplicity, azimuthal 
symmetry in the noise distribution is assumed. 

Since the array beam patterns are not pencil-like with aero response 
in all but the desired direction but are broad and possess a side- 
lobe structvure as shown in Pig. 8, the effect is to smear the fea- 
tui^s of the tru»e field. 

The problem of reconstructing the true from the measured 
distribution is common to the field of astronomy where a radio an- 
tenna is used to survey the distribution of radio brightness over the 
sky. Tlie argments of Bracewell^ ^ who i'irst considered the problem 
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in detail in conjunction with aerial smoothing in radio astronomy may 
be adapted to the underwater acoustics case as follows. Upon inte- 
grating orer the variable, , equation (1) becaaes 

S)X' {^)d& (2) 

where G signifies the average of G over all azimvtth for a given 
colatitude and X = (<^) sirK^*. The integral on the ri^t, a con- 

volution integral, may be written 

Y{s) = Gy^;>^ (5) 

By means of the convolution theorem, equation (5) may be rewritten in 

terms of the Fourier transforms of the three variables as follows 

Y{s) = G{s)X[s) (4) 

where the convolution symbol has been replaced by the algebraic pro- 
duct and the three terms are the spec t rums of the corresp aiding func- 
tions of equation (5) and "s" is the number of cycles per unit of 
angle, ^ . Prom (A) it is seen that K ( s ) = 0 whenever (?(5)=0 
and thus X {s) is indeterminate under such conditions. Obviously 
in any restoration process dependent upon measured values of Y only 
the low frequency components of X are restored to their full value. 
The resulting values of X have been called the "principal solution" 
which may be defined as that solution whose spectrum is the same as 
the spectrum of the true distribution at all values of S for which 
G is] ^ 0, The results of a restoration process to be described in 
the paragraphs following, are shown in Pig. I 5 . The significant 
featxure of Fig. 15 that the transform of the restored field 
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approaches zero at almost the same rate as the transform of the mea~ 
svired field. This is an important condition since any values ofXis) 
for 5 greater than 5 ^ , the value of 5 where G{s)=0 , would rep- 
resent fictitious invisible distributions introduced by the restora- 
tion process. Since « 0 , the addition of fictitious high 

frequency canponents is negligible and X {&) may be considered an 
exact principle solution throughout this study. 

6.2 The Matrix Approximation 

The restoration process ^employed can be> considered in two 
parts! the estimation of a first approximation, and an iterative 
solution deriving from such first approximation. The former, the 
matrix approximation, is novel and is considered in this paragraph. 
The latter, similar to methods of Bracewell, is described in para- 
graph 6.5. 

The matrix approximation is generated by converting the con- 
volution integral of paragraph 6.1 to a finite summation of terms of 
number equal to the number of beams. Each term in the summation rep- 
resents a contribution from a particular solid angle and is propor- 
tional to the average of tlie true field subtended by the solid angle. 
Since a similar expression can be written for each beam, a square 
matrix results relating values of measured field to the unknown 
values of the true field. By inverting the matrix, the true field 
values can be obtained. An exact matrix can be written only for the 
isotropic case. However, the matrix so derived when used with non- 
isotropic distributions yields true field values which although 
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inexact, are closer approximations to the true field distributions 
than can be arrived at by intuition. The approximation so derived 
nay be then conveniently operated upon by iterative techniques to 
yield the principal solution. The details of the derivation of the 
matrix approximation are contained in Appendix III, 

6,^ The Iterative Solution 

The iterative solution is a method of restoring the true 
field by successive substitutions. The first approximation pre- 
viously described is tested by numerically completing the integra- 
tion of cqiaation (1) of paragraph 6,1, and comparing the results with 
the measured values of beam poivor, K, , From experience with nxan- 
erous cemputations of this nature, it has been determined that a 
sviitable second approximation is obtained by correcting the first 
approximation by an amount equal to twice the difference In decibels 
between the measxired and computed values of Y- . Three cuch itera- 
tions were usually sufficient to reduce the difference between mea- 
sured and computed values to less than experimental measurement 
error, 

Bracewell^ has shown that a sufficient condition of con- 
vergence of the above process it that Gks) be non-negative, a condi- 
tion vdiich is satisfied by the transforms of the beam patterns of the 
experimental array, 

A comparison of the decibel differences between a given pair 
of beams before and after restoration show a two to three-fold in- 
crease in the decibel difference. It is reasoned that measurement 



errora will be aimllarly increased In the restoration process. Ac- 
cordingly, the errors in the true field principal solution are es- 
timated to be no greater then * 0,6 decibels except in those regions 
where abrupt discontinuities in the distributions exist. A sample 
computation of the complete restoration process is oontainsd in 
Appendix IV, 

The resxilting true field distributions for the fire environ- 
mental conditions described in paragraph 4.2 are shown and tabulated 
in Figs. 16 through 20, The existence of a marked noise null in the 
horizontal direction is especially ixseful in the design of sonar 
arrays whoso objective is to detect weak si^ials in a noise back- 
groimd. For the detection of horizontally arrtvlng signals in a 
background of noise arriving from all directions, good vertieal 
directivity as well as horizontal directivity is advisable. To 
illustrate} consider two arrays of eqml horizontal directivity but 
differing vertical directivity. Let one array have uniform vortical 
response and the second 20 decibel rejection over all vertical an^s 
except within 15 degrees of the horizontal. The signal-to-noise 
ratio of the latter array will be 5»5 6 decibels better than the 

former when both are positioned in an isotropic noise field and the 
field of Fig. 16 respectively. 

As previously noted, variations in the noise field distri- 
bution for differing array depths could not bo established because 
of uncei^inty as to the constancy of sea state during the recording 
period. The increase in anisotropy with increasing sea state, how- 
ever, is obvious. 




FIGURE 16 TRUE FIELD VERTICAL DISTRIBUTION 
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FIGURE 17 TRUE FIELD VERTICAL DISTRIBUTION 
SEA STATE 3 ARRAY DEPTH 560 FT 
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FIGURE 18 TRUE FIELD VERTICAL DISTRIBUTION 
SEA STATE 3 ARRAY DEPTH 150 FT. 
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FIGURE 19 TRUE FIELD VERTICAL DISTRIBUTION 
SEA STATE I**" ARRAY DEPTH 130 FT. 
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FIGURE 20 TRUE FIELD VERTICAL DISTRIBUTION 
SEA STATE 1/2 ARRAY DEPTH lOOOFT 



VII. AN OCEAN MODEL 



7>1 Ocean Model Deyelopment 

From the nature of the distributions shown in Pigs. 16 
through 20 and oonsistent therewith, the following simple picture of 
the mechaniao of noise generation and redistribution in the ooean 
has been constructed. Noise in the deep ocean is presimed to be a 
superposition of two fields. The first field is that which exists at 
sea state zero prior to the development of wave motion and/or white 
caps. Its distribution is isotropic and its origin was not estab- 
lished during this experiment. A possible explanation of the origin 
of low sea state ambient noise is considered in paragraph 8.1. The 
second source of noise is generated by wave motion and *dilte caps at 
the ocean surface. The equivalent surface radiator is highly direc- 
tional in the vertical plane when azimuthal effects are averaged. 

The directionality can possibly bo attributed to the ocean waves and 
swells functioning in a manner analogous to angular focusing reflec- 
tors. In addition, in azimuth the ocean swells impede transmission 
at shallow transmission angles in directions normal to the wave 
fronts of the swells. 

It is ftirthor postulated that the noise Intensity at a given 
point in the ocean duo to the equivalent surface radiator arrives by 
both a direct ray path and single botcneo path experiencing specular 
reflection at the ocean bottcm. The intensity at the surface source 
experiences inverse square spreading loss and absorptive attenuation 

-of 

of the form lO lO where a is the absorption coefficient 
( decibel s/kiloyard) and R is ray path length in kiloyards. 
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Referring to the model of Fig. 21 eind disregarding for the mcxaent 
the low sea state isotropic contribution to the tote.1 noise field, 
the abore postulation of an ocean model may be restated in mathe- 
matical form as follows: 



X{ff)dn-- 



Si^)d/} 



- aRp 

10 



where X {R)dS2 «= noise intensity per incremental solid angle 

at the receiver in the upper hemisphere 
S{^)dA =■ noise intensity at one yard per incremental 
radiating surface. 



( 1 ) 



Since 



dn = 



d A 



R, 



D 



equation (1) may bo rewritten as 

dS2 lO 



X{R)dS2 






aRp 

lO 



( 2 ) 



( 5 ) 



or, finally, 



q//sec 

S{A) -- X(d) cos ^ lO lO 



0 < (9 < rr / 2 






A similar expraseion for tha intenilty at the receiver in the lowr 
hemisphere due to bottom reflection may be written 






o[Rf{\ + Rff?) 
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lO 



( 5 ) 



where a { d) i.% the intensity bottom reflection coefficient and 
is the supplement 
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FIGURE 21 OCEAN MODEL 
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Since 




(6) 



equation (5) may be written as 





0 < ^ < rr/2. 



(7) 



TT/2. < Spf < V 



1,2 Deduced 5 and a i6) 

Employing the distributions of Figs* 16 through 19 «md equa- 
tions (4) and (7) above, values of 5 (<^) and a (<^)were computed. It 
was necessary to estimate the contribution of the isotropie zero sea 
state noise field and to subtract that amount from the total field, 

X , prior to proceeding with the eomp\itation. Several alterna- 
tive methods of estimating the zero sea state contribution, none of 
»d\leh can be considered exact, are possible. For example, the level 
of the null occurring at an elevation angle of approximately 100 de- 
grees in all distributions may be considered due to the low sea 
state isotropic field. Errors in this assumption are likely duo to 
the difficulty in establishing the true minimum. The most probable 
point of error in the restored distribution is at the minimum due to 
the difficulty Inherent in resolving a sharp discontinuity with a 
broad beam pattern array. In addition, contrary to the very simple 
straight ray assvsnptions of the ocean model, curved path axrivals 
direct from the curface or by convergent zone paths^^^ duo to vorti- 



cal velocity gradients probably contribute to the null level. For 
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example, Bell ^^5) y^as eonaidered the deep ocean ease with negative 
gradient to deep depth. He reports bottcsn reflection loeaee for 
grazing angles of incidence which are approximately 6 decibels less 
due to the refractive effect of the negative gradient than would be 
predicted on the basis of straight ray path assumptions. 

The best estimate of the isotropic contribution utilized the 
Knudsen^^^ sea state/nolse level curves. If the restored distribu- 
tion is averaged over 4wJ? , the relative level derived from an 
omnidirectional hydrophone can bo obtained. Subtracting from this 
level the decibel difference between the Knudsen noise level for the 
sea state under cwisideratlon and sea state zero, the level of the 
isotropic contribution can be established. The weakness in tiie above 
approach rests in its dependence on subjective estimates of sea state 
level. The error introduced by this technique decreases with in- 
creasing sea state since the percentage contribution to the total 
field of zero sea state isotropic noise becomes smaller under such 
oonditimis. 

Figure 22 is a plot of5((^), the svurface source radiating 
pattern, deduced by the above procedures from the restored fields of 
Pigs. 16 through 20, A reasonable approximation to5((^)is the 
funetloni 

90 -(^ 

5(^) = (cos 6^) '0 0 < <? < 70° 

5 (^) = cos^ d 70° < < 90° 



which is shown superimposed on the experimental data 
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FIGURE 22 



SURFACE SOURCE RADIATION PATTERN 



The bottom reflection coefficient (reflected proesure/ 
incident pressure) derived for conditions corresponding to Fig« 22 
are plotted in Fig. 25 together with experimental data and theoreti- 
cal plots from MacKenzie^^^^ and Marsh^^^^, The fit ic considered 
good in view of the Indirect manner in which the data vrere obtained 
and the several assumptions involved in the analysis. The results 
are considered to substantiate the simple ocean model proposed. 
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FIGURE 23 BOTTOM REFLECTION COEFFICIENTS 
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8.1 Conclusions and Disoussion 

Prior to this experiment, as noted in Chapter I, no experi- 
mental data existed to deocribe the directional distribirtion of deep 
ocean ambient noise. Accordingly, stated objectives of this study 
were to investigate anisotropy of ocean ambient noise, to isolate 
vertical and asinxxthal dependencies, to eotablleh the identity of 
sources contributing to non-uniform distributions and to infer the 
radiation patterns of such sources. All the above research objec- 

•a 

tives have been realized by this experiment. Deep ocean ambient 
noise is anisotropic in the frequency band 750 to I5OO cycles per 
second and at sea states above zero. Vertical and azimuthal depen- 
dencies have been Isolated. The surface of the ocean has been 
identified as a major contributor to non-unifom distribxjfbions and 
the radiation patterns of such surface sovirces have been deduced. 

Vertical anisotropy arises from the superposition of a 
noise field due to the radiating ocean surface upon an isotropic 
field present in the absence of surface radiation. The shape of the 
vertical distribution, Figs, I6 throu^ 20, may be compared to the 
form of an arro%»head. The arrow tip is directed normal to the ocean 
surface and represents maximum Intensity arriving from the svirface 
direction. The arrowhead tail indicates noise originating at the 
ocean sxxrface which has been specularly reflected frcan the ocean 
floor. A dip in the arrot;head tall in a direction normal to the 
ocean bottom is evidence of a lox-rer value of bottom reflection co- 
efficient for normally incident sound than for other angles, A 
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narked null at the t.'aiat of the arrowhead, approximately contained 
within 15 degrees of the horizontal, exists in the vertical dis- 
tribution. The magnitude of the null varies between 2 and 12 deci- 
bels belov; the naxinun intensity normal to the surface for sea stat« 
one-half to three respectively. The null is due to tv;o effects. 
First, contributions from the surface source to the horizontal or 
near horizontal distribution arrive from greater distances than 
contributions to the off -horizontal distribirtions and have exper- 
ienced, therefore, greater attenuation from spreading and abaorp- 
tion. Second, the level of radiation from the surface source in 
horizontal directions is more than 20 decibels below the level of 
radiation in the doi-mward directions. 

The non-unifora noise distribution with a niniaua in the 
horizontal direction is of practical importance to the design of 
sonar receiving arrays. The importance of vertical directivity to 
the detection of horizontally arriving signals is apparent. 

The origin of the Isotropic zero sea state noise field v/as 
not established by this experiment. Of the various sources con- 
sidered by Frosch^^^ distant surface motion is considered most 
likely. Although the sea surface nay be flat at the point of mea- 
surement, the noise produced by distant s\urfaoe disturbances can 
propagate very long distances with negligible attenuation due to 
absorption in the 750 I50O cycle per second frequency band ob- 
served. The scattering and reflection processes experienced by the 
sound in reaching the point of observation covild destroy its origi- 
nal character and an isotropio field result. 



Ocean CTrellc are probably responsible for azinubhal oniso- 
tropy by inpedin^ transnission of surface genemied noioe in direc- 
tions nomal to the vaive front of the svrells. 

The intensity at a receiver due to a unit incrcnental sun- 
face radiator of intensity, S' (<$'), at cne yard, can be considered aa 
the sm of intensities due to direct transrdssion and specular bot- 
teta reflection v;ith attenve.tion due to spherical spreading and ab- 
sorption occurring over the separation botvrecn radiator and receiver. 
The quantity, S [ &) ^ is taarkodly directional in the dov.wirard di- 
rection and can be represented approxinately by the expression: 

96-C 

(cos & ) 0 < :? < 70° 

Skd) = , 

cos^^ 70° < 9 < 90° 

The above ocean laodel is not unique. It does, hov;ever, con- 
form to the measured data as evidenced by the reasonable agreement 
between bottom reflection coefficients detomined frora this experi- 
ment and values reported by other investigators. Further additions 
to the model, as for example the consideration of bottonv^surface 
mvdtiple reflections, are obvious. Such refinements would prosxaaa, 
however, a lcnoi-;ledge of bottom topography and array position rela- 
tive to the bottom greater than existed dviring the experiment. 

Although the data upon which the ocean model is based were 
obtained in ocean depths greater than 1C00 fathoras, there is no 
known feature of the model w}dch would preclude its application to 
shallower depths of perhaps 100 fathoras. Its application to oven 



shallower depths will depend upon the presently unknown effect of 
water depth upon the nature of the surface source radiation pattern 
On the basis of neasurementa conducted over a one hour per- 
iod under constant condltlone of array orientation and sea state, 
the ambient noise field can be considered stationary. A single, 
trainable and tlltable dlreotlonal receiver could be employed, 
therefore, to duplicate the performance of the element, multi- 
beam array. However, to achieve beam patterns equivalent to those 
of the 52 element array with a single directional unit, focusing 

-a 

reflectors of dimensions comparable to the J2 element array would 
be required. A very difficult problem of mechanical training and 
elevating the reflecting unit would be substituted for the compara- 
tively simple telemetering problem associated with the 52 element 
array. 



8.2 Suggestions for Further Experiment 

The results contained herein represent only a first step In 
the area of ambient noise directional distribution investigation. 
Objectives for future experiment should include 1 

a. Shallow water investigation. 

b. Extension to higher sea states, 

e. Extension to array depths of at least 5^00 feet. 

d. Extension of bandwidth covered from less than ^lOO 
cycles per second to 5 kilocycles per second. 

e. Investigation of the effect of land masses such as 
islands and the continental shelf and the Influence of surf 
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An additional area of investigation is to substantiate the- 
oretically and experimentally the radiation pattern deduced from this 
experiment of a white cap situated on the crest of a swell, 
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APPENDIX 1 



QUANTITATIVE SYSTEM AI7ALYSIS 

The objective of this aeotion is to derive an expreseion for 
the voltage appearing at the output of a given beam detector for a 
rand Ota noise inpxit. 

Consider first an N element array of omnidirectional elements 
situated in an isotropic noise field of purely continuous spectrum. 

If the hydrophone electrical outputs generated by the noise field are 
summedi with or without the introduction of an electrical delay, and 
if the hydrophones are sufficiently widely spaced such that cross- 

correlations between elements are zero, the array output following 

1/2 

square law detection will be proportional to .If the cross- 

correlations are not equal to zero, the otitput will be proportional 

to A/ ( I + 6’ ) ^ where 

I N^N-W 

^ = 77- (0) 

and /o is the normalized cross-correlation function derived from 

</ 

the input noise spectrum. The value of e may also be determined 
from the directivity ratio. Utilizing the beam patterns of paragraph 
2.6, the directivity ratio is derived by nmerlcal integration from 
the expression 




where G (<^ij^)ls the beam pattern of paragraph 2.6, and ClOjO)- I . 
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The quantity, " e" , may be determined from the expreasion 
/V / \ + e - \ / D. R t For inner and outer beams " s" has been ccxn- 
puted to be O.25 and 0,02 respectively, Lawson^^^^ hag shown that 
the correlation functions of imolipped and clipped signals are re- 
lated by p -2./v arcs in P * snail values of correlation 

■xlipped 'J unclipped 

as is the case under consideration, the arcsine of the argvsnent is 
equal to the argvonent and the correlation function of the clipped 
signal reduces to 2 / tt times the imclipped correlation function. 
Since the actual detector employed was a linear rectifier rather than 

■a 

a square law device, an added factor K > 0,8 is required, K -'J2/W 
is the ratio of the mean absolxxte value to the root-moan- square 
value in a normal distribution, A further factor of 0,5 ie intro- 
duced because the scale of the abscissa of Fig, 11 is I 6 for maximxa# 
detector output rather than J2, The detector output voltage is 
therefore proportional to 




Substitxrtion of the values A/ ■ J2 and e ■ 0,25 into the above ex- 
pressions produces the value 2,44, This quantity when entered as 
the abscissa of Fig, 11 indicates an output value of 52 millivolts 
for detector No, 1, an inner beam. 

The computed value of 52 millivolts was experimentally con- 
firmed from data recorded at sea. When the outputs of all J2 beams 
are averaged, the resulting value is eqxilvalent to an individual 
beam output to be expected fr(xa an isotropic, random noise field. 

The averages of several sets of data eqxialled approximately 52 mlllt* 



volts 



APPENDIX II 



SYSTEM LINEARITY 



Rudniek^^^^ includes the 



necessary expressions to compute 



beam output power as a fixnctlon of array geometry and input slgnal- 
to-noise ratio at an individual hydrophone. Output power may be 
expressed as 



with the summations to be completed over the N array elements. The 
quantity, ( 0 ) , is the combined autocorrelation fvmction of 

signal and noise. For purpose# of this analysis, the spectra of 
algnal->to-nolso may bo considered identical and are assigned dif- 
ferent nomenclature only to differentiate between that noise which is 
within the main lobe of the directed beam and that which is outside 
the lobe. Implicit in the last statement is the assumption that 
noiee is completely uncorrelated with respect to angle of arrival, 

C !j ( 0 ) is the autocorrelation function of the spectnan introduced 
by clipping and is defined by Rudnlck^^^^ as 




( 1 ) 




( 2 ) 



ju. ( 0 ) may be expanded to 
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where p and a are noise and signal autocorrelation functions re~ 
spectlrely, (O) = I since by definition it is the signal to 

-V V 

which the beam is steered, S C^,(0)mB.y also bo separated into 

/ti / --I 



/'/ 
/- 1 



A/ N N[N-\) 

Z C -- I cAO)^ II 

j-\ /-I 



' 



( 0 ) 



(A) 



It follows from (2) that = ^/2 -| • In addition, Rudnick has 

shown that 

AJ N 

II c A(S) 

is small and may be neglected. 

Substituting the several expressiona described above into 
(1) and reducing terms, the output voltage of the linear half-wave 
rectifier can be shown to be proportional to 




2 

7T 



n 

n->rS 





(5) 



If a small signal-to-noise ratio input of s/n equal to 0,1 
is assumed and with N equal to 52 and "e*" equal to 0,25, equation 
(5) yields the value ^,61. Entering the detector characteristic of 
Fig, 11 with the above value, an output voltage of 90 millivolts is 
obtained. No recorded data was used in the analysis where the out- 
put of any beam exceeded 90 millivolts. It is concluded, therefore, 
that the system was operated throughout its linear region diorlng the 
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APPENDIX III 



DERIVATION OF THE MATRIX APPROXIMATION 



Considering the case of aziauthal symmetry, the oxrtput of 
any Tertleal beam may be written* 

y, -- J X{^) G, dn (1) 

•^477 

where * output of / beam 

X » continuous true field distribution (power/unit 

receiving area/untt solid angle subtended by the 
source) 

G “ power response of / beam 
If the true field is assvmied to be isotropic, then 




The directivity ratio of an array is defined as* 



( 2 ) 



G, {^,<^)dI2 

7T 

G, {^, f dn 

Combining equations (2) and (5) 




(5) 



Y, = X^ G %cp) D.R ^7T (4) 

where G equals unity and 0.855 for outer and inner beams 

respectively as defined by paragraph 2.6. Setting 4 tt 
which has the physical Interpretation of total power from all direc- 
tions normalised to unity, then 
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D.R, (outer beam) 

>: = ( 5 ) 

0,855 D.R. (inner besim) 

Rewriting equation (I) ae a sum of eight terms, with each term rep- 
resenting the contribution from a spherical aone formed by two 
planes parallel to the equator and each zone eontaining one of the 
eight vertical beams, thenj 

■' = 8 . Zrr 

X. I sin (Pdi9 

^ •'zone 

The true field, X j , may be non-isotropic, although X j is a constant 
within a particular zone. The quantity G ^ may next be aver- 

aged over all azimuth. 



J G^ { & , <^ ) d(p 



( 6 ) 



Jz\ 



K = I' K 2^ G s\x\ddd 



( 7 ) 



where d and & are the limits of the spherical zones. There 
exists between the limits of each zone a particular value of<^=<^,p, 
which permits equation (7) to be transformed to the following 



/ = 8 



Y, --Z 

J- I 





( 8 ) 



Equation ( 8 ) may bo further simplified by replacing the integral by 

a constant, k ^ , for a particular zone yielding 

J-B 



Y. ^ 2 ttG, { ) Z 

J ~ ^ 



(9) 
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or more briefly 

y, = Xj /, y = 1-8 (10) 

Equation (10) relates the measured field Talues^ , to a sot of 
eight simultaneous equations in eight \xnknowns vhich are the true 
field Talues, Xj . 

The a ,j matrix which satisfies equaticxis ( 5 ) and ( 9 ) was 
determined by trial and error. The matrix was then inverted by 
digital computer to produce a set of equations of the form 

The coefficients of the inverted matrix are listed below. By mul- 
tiplying each element of row one by , row two by Y^ I etc , p the 
sum of all column one elements yields X^ , the column two elements 
X^, etc. The values ot X so derived are the first approximation 
to the true field in the restoration process. 

Table of Coefficients < 7 ^^ 



V 


1 


2 


5 


4 


5 


6 


7 


8 


1 


67.855 


-52.557 


27.269 


-19.912 


10.678 


-19.522 


5.122 


-8.952 


2 


- 1 . 67 ? 


89.527 


- 67.679 


25.059 


-19.5^ 


16.75^ 


-5.158 


- 1 .845 


5 


-0.056 


- 9.622 


92.556 


-59.901 


20.905 


-52.688 


0.808 


- 0.282 


4 


- 1.102 


2.505 


-52.047 


106.681 


- 58.812 


20,422 


- 5.628 


- 0.799 


5 


- 0.799 


- 5.628 


20.422 


-58.812 


106.681 


-52.047 


2.505 


- 1.102 


6 


- 0.282 


0,808 


-52.688 


20.905 


-59.901 


92.556 


- 9.622 


-0,056 


7 


- 1 .845 


-5.158 


16.75^ 


-19.5^5 


25.059 


-67.679 


89.527 


-1 .677 


8 


-8.952 


5.122 


-19.522 


10.678 


-19.912 


27.269 


-52.557 


67.855 




( 






APPE>nDIX IV 



TRl’E FIELD RESTORATION SAXPLE COMPUTATION 

In order to Illustrate the restoration process, the proce- 
dures followed In deriving Fig. 18 will be described. 

Table I ie a random selection of recorded data typical of 
that which results following the film reading procedures of para- 
graph 5.6. The values shown are the Indlvldml beam detector out- 
puts in millivolts following one second Integration and at 14 second 
Intervals. Beam numbering is such that beam is directed at the 
surface, beams 27-J1 at J8 degrees from the vertical, beams 22-26 
at 64 degrees from the vertical, beams 17~21 at 78»5 degrees from 
the vertical, etc., to beam 1 directed at the bottom of the ocean. 

The RMS values for each beam for each rim are computed and tabulated 
in colvsati 2 (Aver) of Table II, The detector correction is applied 
(Corr Aver) and the RMS output for beams of a given colatitude com- 
puted (Vert Aver), The average of all beams is similarly determined. 
Power ratios for each colatitude with respect to beam 52 and 
with respect to the average of all beams are computed. The Y, 
column is determined by multiplying the appropriate expression of 
equation (5) of Appendix III by the power ratios of column 8 of 
Table II, The values of Y, are then used in conjunction with the 
matrix of Appendix III to obtain the first approximation to the , 
In Fig, 24 the measured values of , expressed as decibel 
ratios with respect to beam 52, are plotted as dots, A faired 
curve connecting the first approximation values of is indicated 
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Beam 

1 

2 

5 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

15 

14 

15 

16 

17 

18 

19 

20 

21 

22 

25 
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TABLE I 

51 50 51 50 51 51 51 52 50 51 51 50 50 50 51 50 51 51 50 

49 49 49 48 48 48 49 49 48 49 48 48 49 49 48 4<5 47 49 49 

47 46 4€ 46 47 46 47 48 47 45 46 47 47 46 4^ 4^ 46 47 47 

49 49 48 48 49 48 48 49 49 48 49 49 49 49 49 48 49 49 49 

50 50 49 49 50 49 50 49 49 49 49 49 50 50 49 49 50 51 49 

47 47 47 47 47 46 47 46 45 47 47 46 46 47 46 46 47 46 47 

49 50 49 49 49 48 49 50 49 54 51 50 50 49 50 49 50 50 50 

50 51 50 50 49 50 50 4l 50 50 50 50 51 ^9 50 50 51 50 50 

48 49 49 48 49 49 50 49 48 49 49 49 49 48 48 47 4o 50 48 

49 52 50 50 50 50 51 52 50 50 50 51 50 50 50 51 51 52 50 

46 48 47 46 47 46 48 47 48 46 47 47 48 48 47 47 48 48 46 

42 44 45 44 45 44 45 45 45 4j 44 44 45 44 44 45 44 45 45 

42 44 45 45 45 42 44 45 44 45 44 45 4j 4j 4j 45 44 45 44 

4l 45 45 45 45 42 44 45 45 4l 42 4j 42 42 42 42 4j 44 4j 

46 49 49 49 49 48 49 50 48 47 49 49 49 49 48 47 48 49 47 

50 55 54 54 54 51 55 54 55 50 5^ 51 52 52 52 52 51 52 52 

47 50 49 49 49 49 49 49 48 56 52 50 50 49 50 49 49 50 49 

54 58 51 56 56 54 57 5^ 56 56 57 5^ 5^ 55 55 5^ 5<5 5^ 5^ 

48 51 51 50 50 50 51 50 50 49 49 50 51 50 49 48 50 52 49 

50 52 51 50 51 50 51 5^ 52 51 50 50 49 50 51 51 50 50 51 

54 56 51 55 57 56 55 55 55 52 5^ 5^ 55 55 55 52 55 5^ 5i 

50 54 51 51 51 50 52 55 54 52 52 55 52 52 52 50 52 52 51 

64 64 64 65 64 65 65 65 64 61 66 65 66 64 60 65 65 64 65 



Beam 

24 

25 

26 

27 

28 

29 

50 

51 

52 
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TABLE I (CONT.) 



62 62 62 61 61 61 62 61 60 59 61 62 62 60 62 59 60 
51 55 50 51 50 50 50 52 51 50 52 51 52 51 51 51 51 
57 59 58 58 5<5 5^ 5<5 59 57 5<5 58 57 57 58 58 57 57 
60 59 59 60 58 58 59 58 58 56 58 58 59 59 58 57 58 
60 59 60 61 60 59 59 59 61 58 60 61 60 58 58 57 53 
66 62 62 65 62 64 65 64 65 62 65 66 62 62 64 60 6i 
69 74 67 65 64 64 66 65 61 64 66 69 67 68 65 64 70 
64 65 62 64 61 64 62 66 62 60 62 64 61 65 62 62 61 
77 78 76 74 74 76 75 77 75 74 75 76 72 74 75 76 77 



65 59 
51 51 
60 57 

60 57 
62 58 
65 65 
65 65 

61 59 
75 72 
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TABLE II 







Corr 


Vert 


Elev 


Ratio 


Ra^tio 


Ratio 


Ratio 




Bead 


Aver 


Aver 


Aver* 


Angle 


re 52 


re av. 


Re 52 (db) 


Re av.(db) 




1 


52 


52 


52 


180 ® 


.605 


.999 


-2.2 


0.0 


.0555 


2 


47 


49 












5 


46 


48 
















4 


49 


51 


50 


142 ® 


.549 


.909 


-2.61 


- 0.42 


.0287 


5 


48 


52 














6 


45 


48 
















7 


4 i 


4 o 
















8 


55 


50 
















9 


48 


47 


47 


116® 


.495 


.821 


-5.06 


-0.87 


.0276 


10 


48 


51 












11 


45 


47 
















12 


45 


45 
















15 


4 i 


55 
















14 


44 


47 


46 


101.5® 


.475 


.784 


- 5.26 


-1 .06 


.0248 


15 


47 


49 














^6 


44 


44 
















17 


55 


55 
















18 


47 


48 
















19 


45 


47 


52 


78 . 5 ® 


.611 


1.010 


-2.15 


0.05 


.0519 


20 


46 


48 










21 


61 


62 
















22 




56 
















25 


5 ^ 


52 
















24 


52 


52 


55 


64 ® 


.671 


1.122 


-1.74 


0.46 


.0576 


25 


55 


— 
















26 


61 


59 
















27 


54 


58 
















28 


54 


58 
















29 


— 


— 


60 


58 ° 


.797 


1.520 


- 0.98 


1.21 


.0417 


50 


59 


64 
















51 


57 


59 
















52 


^7 


67 


67 


0 ° 


1.000 


1.655 


0 


2.19 


.0556 



•Vertical Average - The averaged output of all beams of eqml eleva- 
tion angle. 
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FIGURE 24 TRUE FIELD VERTICAL DISTRIBUTION 
HERATIVE SOLUTION 



by dashed lines. This curve is used to select values of Xj in 10 
degrees steps to complete the numerioal integration of equation (1) 
of paragraph <5.1, The integration of equation (l) was accoaplished 
with incremental steps of 10 degrees in elevation and aulnuth. 
Figxxro 25 represents the results of all eonputations with the ex- 
ception of the integration over the variable, d . The ordinate is 
G^d) Xi6^)s\n^, the abscissa, 9 . The 8 curves represent the 8 
vertical beams and is included to illustrate the fact that the beam 

response contributes significantly in all directions and emphasizes 

% 

the importance to this process of knowing the bean patterns exactly 
The curves of Fig. 25 are then integrated by planimeter sind the 
decibel ratios between a given beam and the beam directed at the 
surface obtained. The results are plotted as «Vsin Fig. 24. The 
decibel difference between dots and X's is moasxired and twice that 
value applied as a correction to the first assumption, dashed curve 
of Fig. 24. The dotted curve resvdts. The numerical integration 
process is repeated yielding the open circles. A final iterati(xi 
yields the solid curve which is the curve of Fig. 18. 




FIGURE 25. G{e)X{e)S\N$ VS B 
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